Int. J. Mol. Set 2012, 13, 5584-5597; doi:10.3390/ijms 13055584 



OPEN ACCESS 



International Journal of 

Molecular Sciences 

ISSN 1422-0067 

www.mdpi.com/journal/ijms 

Article 

Loss of Genetic Diversity in the Cultured Stocks of the 
Large Yellow Croaker, Larimichthys crocea, 
Revealed by Microsatellites 

Le Wang , Xiaofeng Shi 2 , Yongquan Su 2 , Zining Meng '* and Haoran Lin *'* 

1 State Key Laboratory of Biocontrol, Institute of Aquatic Economic Animals and the Guangdong 
Province Key Laboratory for Aquatic Economic Animals, School of Life Sciences, 
Sun Yat-sen University, Guangzhou 510275, China; E-Mail: lewang.wang@hotmail.com 
College of Ocean & Earth Sciences, Xiamen University, Xiamen 361005, China; 
E-Mails: xiaofengshi@xmu.edu.cn (X.S.); yqsu@xmu.edu.cn (Y.S.) 

* Authors to whom correspondence should be addressed; E-Mails: mengzn@mail.sysu.edu.cn (Z.M.); 
lsslhr@mail.sysu.edu.cn (H.L.); Tel.: +86-20-84110188 (Z.M.); +86-20-84113791 (H.L.); 
Fax: +86-20-84113717 (Z.M.; H.L.). 

Received: 14 March 2012; in revised form: 23 April 2012 / Accepted: 2 May 2012 / 
Published: 9 May 2012 



Abstract: The large yellow croaker (Larimichthys crocea) is the most important 
mariculture fish species in China and the wild stocks of this croaker have collapsed in the 
past decades due to high fishing pressure and habitat degradation. Due to a lack of wild 
croaker samples, however, studies concerning the genetic changes of the cultured croaker 
stocks compared to their wild counterparts were never conducted. Here, we collected three 
wild populations in the northern and central East China Sea during fisheries survey and 
investigated the differences in terms of genetic diversity and differentiation between and 
within cultured stocks and wild populations. Our results demonstrated that the cultured 
croaker had significantly reduced genetic diversity in contrast to the wild populations, and 
also presented statistically significant differentiation from the wild, indicating that 
enhancement of the current wild stock should be conducted with caution. These changes 
may be caused by founder effects, artificial selection and random genetic drift. With a 
relatively high level of genetic diversity, the wild populations showed important value for 
improving the ongoing breeding program of this croaker. Further, we detected no 
differentiation among the wild populations, suggesting that the wild croaker in the northern 
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and central East China Sea should be considered as one unit for management 
and conservation. 

Keywords: large yellow croaker; Larimichthys crocea; genetic diversity; differentiation; 
micro satellites 



1. Introduction 

The large yellow croaker, Larimichthys crocea, is an economically important marine fish species 
endemic to China. The production of this croaker reached about 200,000 tons in the mid 1970s and it 
was once ranked in the top three commercial fish species in mainland China [1,2]. However, the wild 
resources of this croaker has collapsed in the past decades due to heavy exploitation of spawning and 
over- wintering aggregations, poor stock management, habitat pollution and climate changes [2-4]. 
In order to satisfy the needs of consumers for food and also to protect this species from extinction, the 
Chinese government has conducted successful artificial mariculture for this croaker since 1985 [5]. 
The aquaculture production reached approximately 70,000 tons in 2006 [6]. 

However, aquaculture practices are likely to reduce the genetic diversity and further to cause the 
loss of disease resistance and environmental adaptability, which greatly limit the potential for selective 
breeding [7-9]. Recently, several biological changes including small size and early age of sexual 
maturation, low growth rates, poor flesh quality and loss of resistance to disease and cold have been 
identified in the cultured croaker in contrast to the wild populations [2,6]. Such changes were 
suggested to be associated with overexploitation of the wild stocks and mariculture operations [2,6]. 
The decline of quality may be caused by loss of genetic diversity in the cultured stocks [8]. Following 
successful hatchery production, larvae and fingerlings from hatcheries were released by the Chinese 
government to restore and enhance the wild stocks of this croaker [10-12]. Such artificial release and 
random escape from mariculture stocks into the open marine environment can cause potential harmful 
effects on the genetic make-up of both wild and reared populations without screening the genetic 
backgrounds of the two types of stocks [7,13,14]. However, no studies have been performed to monitor 
such genetic make-up and changes between cultured stocks and wild populations of this croaker 
because of great difficulties in collecting enough wild samples for population genetic studies. 

We successfully collected wild populations of the large yellow croaker in fisheries survey from 
2007 to 2010, which allowed us to conduct genetic studies as described above. In addition, microsatellite 
markers have been successfully used in genetic monitoring of the changes between hatchery stocks and 
wild populations [15-17]. Here, we used 10 microsatellites to analyze the genetic status of both 
cultured and wild populations of the large yellow croaker in China. The aim of our study was to 
examine the changes of genetic variation and to assess potential genetic differentiation between 
cultured and wild populations of this croaker. 
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2. Results 

2.1. Genetic Variation within Populations 

Under exact tests, no consistent deviations from HWE in each sample were detected. After 
sequential Bonferroni correction, only two locus-population pairs (H37-NB and H54-ZJW) were 
significant (corrected P = 0.0007; Table 1). We also found no evidence of LD among loci. Significant 
presence of null alleles was detected by Microchecker at one locus, H43. The estimated null allele 
frequencies in each sample at this locus were more than 0.1 and were statistically significant. As the 
presence of null alleles may bias the results of genetic variation and particularly differentiation, we 
therefore excluded this locus from further analysis. 

Table 1. Summary statistics at nine microsatellite loci across samples of large yellow croaker. 
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Table 1. Cont. 
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* Significant deviation from Hardy- Weinberg equilibrium (HWE) after Bonferroni correction (P < 0.0007); NA, not available. 

In terms of genetic variation, one locus, H80, was monomorphic in two cultured stocks, DQ and 
MY (Table 1). The other loci all showed a high level of polymorphism in each sample. The highest A 
was observed in the wild population YCW (7.899), while the lowest was in cultured stock DQ (4.444). 
For A R , the highest and lowest values were also found in wild population YCW (7.816) and in cultured 
stock DQ (4.307), respectively. In terms of heterozygosities including Ho and H E , the highest and the 
lowest values were also identified in wild population and cultured stock, respectively (Table 1). 
Further investigation revealed a statistically significant reduction of genetic diversity in the cultured 
stocks compared to the wild populations. In detail, A R varied from 4.307 (DQ) to 5.037 (NB) with a 
mean of 4.686 in cultured stocks, which was significantly lower than in wild populations ranging from 
5.913 (ZJW) to 7.816 (YCW) with a mean of 6.730 (P < 0.01). Similar to Ar, H 0 and H E values were 
also significantly lower in the cultured stocks (from 0.427 to 0.491 and from 0.462 to 0.517 for Ho and 
H E , respectively) than in the wild populations (from 0.570 to 0.647 and from 0.591 to 0.649 for Ho and 
H E , respectively, P < 0.05). However, Fis values showed no significant differences between the 
cultured stocks and the wild populations (P > 0.05). The Bottleneck analysis did not detect signals of 
recent population reduction for each sample whether in test under TPM or in mode- shift test. 

2.2. Genetic Differentiation among Populations 

The global Fst was 0.070 with 95% confidence interval ranging from 0.023 to 0.173 and was highly 
significant (P < 0.001). Pairwise Fst analysis showed that the wild populations were significantly 
divergent from the cultured stocks after Bonferroni correction for multiple comparisons (Table 2; 
P < 0.002). Among five cultured stocks, only DQ was significantly different from the others. 
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Interestingly, there was no differentiation detected among the three wild populations (Table 2). Results 
of AMOVA further supported the significant genetic differentiation between wild populations and 
cultured stocks, which occupied 10.208% of total genetic variation (P < 0.001; Table 3). FCA analysis 
also revealed significant differences in terms of allele frequency between the cultured stocks and 
the wild populations (Figure 1). These results were in accordance with the pairwise Fst analysis. In 
simulations of the Bayesian approach with the program Structure, the mean Ln Likelihood values 
clearly suggested three clusters as the most likely population structure (Figure 2). The results indicated 
that almost all the wild individuals were assigned into one cluster, whereas the cultured stocks showed 
identical genetic properties (Figure 3). However, we observed that the cultured stocks and wild 
populations were clearly assigned into their own clusters at K = 2 (Figure 3). These results strongly 
supported the results of the other genetic differentiation studies. In addition, we also observed that the 
DQ stock was slightly divergent from the other cultured stocks both in FCA and in Structure analysis 
(Figures 1 and 3), which was consistent with F S t results (Table 2). 

Table 2. Pairwise F S t (below diagonal) and associated P values (above diagonal) among 
cultured stocks and wild populations of the large yellow croaker. Significant P values after 
Bonferroni correction for multiple comparisons are denoted in bold. 
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Table 3. Results of analysis of molecular variance (AMOVA) between cultured stocks and 
wild populations of the large yellow croaker. 



Source of Variation Sum of Squares Variance Components Percentage Variation P Value 
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Figure 1. Three dimensional scatter plots for individuals of eight samples of the large 
yellow croaker based on factor correspondence analysis (FCA) analysis, in which each axis 
represents one principal factor. 




Figure 2. The log likelihood over 20 runs for each K values, where the highest value is 
suggested to be the true number of clusters (K = 3). 
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Figure 3. Results of structure analysis based on nine microsatellite loci. Each individual is 
represented by a vertical line, which is colored according to the assigned groups at 
estimated K = 3 (see Figure 2). Results of K = 2 are also presented for identifying the 
significant differentiation between cultures stocks and wild populations. 
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3. Discussion 



3.1. Micro satellites Polymorphism 

Investigation of genetic variation of cultured and wild populations of domesticated animals can 
provide valuable information for breeding programs and also for conservation genetics [8,18]. In this 
study, we used polymorphic microsatellites to analyze the genetic differences between the cultured 
stocks and the wild populations of the large yellow croaker. The allele numbers per locus varied from 
3.125 to 16.125 and from 0.000 to 0.889, respectively, which is similar to other marine fish species, 
such as Atlantic cod (Gadus morhua), sea bream (Pagrus major) and the orange-spotted grouper 
(Epinephelus coioides), suggesting these microsatellites are sufficient to detect genetic variation in the 
large yellow croaker [19-21]. Across the eight samples, however, the average allele number and He 
per locus were 5.750 and 0.541, respectively, which is much lower than that (20.6 and 0.79 per locus) 
found in marine fishes [22]. This result might indicate that the large yellow croaker has reduced in 
genetic diversity due to high fishing pressure and/or artificial breeding. 

3.2. Genetic Variation within Populations 

Aquaculture practices have been broadly reported to have the tendency to reduce genetic variability 
in cultured stocks of fish species [16,17,23]. Aquaculture practices are detrimental to the domestication 
process of cultured stocks because a high level of genetic variation is related to adaptive fitness and 
therefore can provide more chances for organisms to survive under the pressure of artificial and natural 
selection [24]. In China, the cultured individuals of the large yellow croaker were consistently 
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considered to show a picture of small size, early age of sexual maturation, low growth rates, poor flesh 
quality and loss of resistance to disease and cold, which may be a reflection of genetic diversity 
loss [2,6]. As expected, in our study, we detected a statistically significant decline of genetic variation 
in the cultured stocks compared to their wild counterparts of this croaker (Table 1). The genetic 
variation measures used in our study involved A R , Ho and H E . Typically, A R is independent of sample 
size and is more sensitive to be detected in populations of reduced genetic diversity than 
heterozygosity, as loss of rare alleles shows little effect on heterozygosity [15,18,25]. Here, we found 
genetic variation declining not only in terms of A R but also in terms of heterozygosity measures, which 
likely demonstrates that the reduction of genetic variation in the cultured stocks is not caused by sampling 
and genotyping bias. 

The decline of genetic diversity in cultured stocks is typically considered to be the result of 
interacting founder effects, random genetic drift, and artificial and natural selection in the cultured 
environments [8,26,27]. In our study, the reduced genetic variation in croaker stocks could be mainly 
due to founder effects during the domestication process. The broodstocks of this croaker were initially 
set up in the mid-1980s, when the wild stock of this species had collapsed [2,6]. In this situation, the 
broodstocks were mainly from random capture of the wild individuals and therefore consisted of only a 
small number of individuals. This type of broodstock is prone to have great effects on genetic variation 
of their offspring and ultimately lead to the loss of genetic diversity in the cultured stocks [27]. 
In addition, random genetic drift and artificial selection much likely played important roles in leading 
to the loss of genetic diversity in the cultured stocks. On one hand, as an endangered species, there is 
no large source population available to supplement the small-sized broodstocks of this croaker. In this 
case, genetic drift is seldom avoided and the cultured stocks would experience excessive loss 
of genetic variability [8]. On the other hand, these five cultured stocks have experienced artificial 
breeding for two to three generations. During such domestication process, artificial selection cannot be 
avoided in order to obtain fingerlings of high quality. Artificial selection has been broadly reported to 
possibly reduce genetic diversity [18,26,28]. Apart from these factors, natural selection and mutation 
may also have important effects on genetic diversity [8], though we cannot clearly differentiate 
whether the two factors reduce or increase the level of genetic diversity in our work. However, it 
should be noted that all the cultured samples have not experienced bottleneck. This may be due to the 
short domestication history of this croaker. 

3.3. Genetic Divergence among Populations 

With respect to genetic divergence, pairwise F S t, AMOVA, FCA and Structure analysis 
consistently supported the significant differentiation between cultured stocks and wild populations 
(Tables 2 and 3, Figures 1 and 3). Such significant differentiation between cultured stocks and wild 
populations was also observed in many other food fish species, such as salmon (Salmo salar), grass 
carp (Ctenopharyngodon idella) and orange- spotted grouper, which was considered to result from 
artificial selection, founder effects and genetic drift [15,21,29]. Interestingly, we did not detect the 
signals of differentiation among the three wild populations (Table 2, Figures 1 and 3), suggesting high 
gene flow in this croaker. It is common for marine fishes, especially for migratory species 
to present little divergence because of lacking clear barriers to gene flow within open marine 
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environments [30]. This result provides important information for wild stock management and 
conservation of this croaker. Among the cultured stocks, we observed that only DQ stock was 
significantly divergent while the others showed little differentiation from each other (Table 2, 
Figures 1 and 3). Considering the fact that all stocks except for NB experienced selection, the shallow 
differentiation among cultured stocks might indicate the effects of artificial selection have not been 
high enough to be detected. 

3.4. Implications for Artificial Breeding and Conservation 

In total, we detected statistically significantly less genetic diversity in the cultured stocks of the 
large yellow croaker in contrast to their wild counterparts, which was mainly caused by founder effects, 
random genetic drift and artificial selection involved in aquaculture practices. In addition, significant 
genetic divergence between cultured stocks and wild populations was also observed. Possessing a high 
level of genetic variation, the wild croaker showed great value for ongoing selective breeding 
programs by providing more genetic variation. At the same time, due to the lack of differentiation, the 
wild populations can be considered as one unit for conservation. However, ongoing artificial wild 
stock enhancement by releasing of cultured croaker fingerlings should be conducted with caution, as 
previous studies have suggested that release of hatchery fingerlings with low genetic variation would 
likely reduce the genetic diversity of the wild populations and further lead to loss of adaptation to 
variable environments for the wild populations [31-33]. The stock enhancement programs of this 
croaker were mainly carried out in the north coast of Zhejiang province and more than ten million 
cultured fries were released from 2000 to 2009, among which about 60,000 were tagged by hanging 
scutcheon [34]. Recapture study suggested that the released croaker could not only survive but could 
also spawn in the following years, although the survival rate was quite low [11,34]. Therefore, 
it is very likely for the released individuals to be integrated into the wild croaker populations. In our 
study, we detected that the cultured croaker was less diverse in genetic variability than the wild croaker 
and was significantly divergent from the wild. If the stocked croaker is incorporated into the wild 
populations, it would definitely change the genetic make-up of the wild croaker populations and 
further cause adverse effects on the wild croaker. 

Since the 1990s, sample collections of the wild populations of this croaker have rarely been 
reported. In our study, we collected three samples of the wild croaker in the northern and central 
East China Sea during fisheries survey for several years, which may suggest that this area is "refugia" 
of this species. Although the wild stock of this croaker has collapsed, we did not detect the signals of 
recent bottleneck in the three wild croaker populations using microsatellites. In another study, we 
analyzed the phylogeography of this croaker using the same wild samples as in this study by 
sequencing mitochondrial Cytb and COI genes. The results revealed that this species was in the 
process of population expansion after Pleistocene glaciations, which also suggested few signals of 
genetic bottleneck in the wild croaker populations [35]. However, previous study has demonstrated 
that the potential of detecting bottleneck is greatly limited in populations experiencing expansion [36]. 
Combining the fact that wild croaker was seldom captured, the effective population size of the wild 
croaker can be rather small. Thus, strict measures by the government, such as reducing fishing pressure 
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and avoiding marine environment pollution in this area must be taken into account to protect this 
species from extinction. 

4. Experimental Section 

4.1. Sample Collection 

We analyzed five cultured stocks and three wild populations of the large yellow croaker. Among 
these five cultured stocks, three (DQ, MY and NB) were collected from Ningbo, Zhejiang province, 
while the other two (XP and SD) were from Ningde, Fujian province. In detail, two (MY and NB) and 
three stocks (DQ, SD and XP) were the third and second generation offspring of the local fish farms, 
respectively. The broodstocks of these cultured stocks were founded using brooders collected from the 
central East China Sea. Due to stock collapse in the wild resources of this croaker, however, it is very 
difficult to capture the wild fish. Fortunately, we collected three wild populations (ZJW, QDW and 
YCW) from the northern and central East China Sea during several fisheries surveys between 2007 and 
2010. Judging from the body size, all the wild individuals were adult. Detailed information about 
sample size and location is shown in Figure 4. A small piece of muscle tissue or fin clip was collected 
and stored in 95% ethanol for DNA extraction. 

Figure 4. Maps for sampling localities of the large yellow croaker used in this study. YCW: 
wild population from the northern East China Sea; QDW and ZJW: wild populations from 
the central East China Sea. DQ, MY, NB: cultured stocks collected from Ningbo, Zhejiang 
province; XP and SD: cultured stocks from Ningde, Fujian province. Sample size is 
showed in number after name abbreviation. 
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4.2. Molecular Methods 

Genomic DNA was isolated using the standard phenol-chloroform extraction protocol [37]. All 
sampled individuals were genotyped using ten micro satellite loci, namely, H16, H31, H33, H37, H43, 
H47, H54, H65, H80 and H82 [38]. For micro satellite genotyping, forward primers were 5'-labeled 
with a fluorescent dye HEX or 6-FAM. PCR amplification was performed according to the Molecular 
Ecology Resources Primer Development Consortium [38]. PCR products were separated on an ABI 
PRISM 3730 DNA automated sequencer (Applied Biosystems) and were measured according to the 
ROX-500 standard using GeneMapper (Applied Biosystems). 

We employed the number of alleles (A), observed heterozygosities (Ho), expected heterozygosities 
(He) and allele riches (Ar) to measure the genetic variation in each sample. These parameters were 
calculated in FSTAT version 2.9.3.2 [39]. Exact tests of Hardy-Weinberg equilibrium (HWE) in each 
sample for each locus and linkage disequilibrium (LD) between pairs of loci were tested using the 
Markov chain methods implemented in Genepop 4.0 [40]. Inbreeding coefficient (F IS ) was also 
estimated using FSTAT version 2.9.3.2 [39]. The genotyping errors and presence of null alleles were 
checked by using the program Microchecker [41]. The significance of differences in genetic variation 
between cultured stocks and wild populations were tested using Mann-Whitney U test. Recent 
population size reduction was examined in the form of heterozygote excess using the program 
Bottleneck version 1.2.02 [42]. The possibility of recent bottlenecks was tested under two-phase model 
(TPM, with 90% stepwise-mutation) with 1000 iterations and using the graphical mode-shift test by 
Luikart et al. [43]. 

Genetic differentiation among samples was estimated using pairwise Wright's F-statistics (Fst)- 
The significance was tested by a permutation with 10,000 replicates using ARLEQUIN 3.5 [44] with 
sequential Bonferroni correction at the significance level of 0.05. Analysis of molecular variance 
(AMOVA) was performed to partition genetic variance hierarchically between the wild populations 
and cultured stocks using ARLEQUIN 3.5 [44]. We also performed factor correspondence analysis 
(FCA) using the program GENETIX 4.05 [45] to detect population structure based on allele 
frequencies. In addition, a Bayesian method was employed to investigate the population structure of all 
samples using the program Structure 2.2.3 [46]. This program can estimate the number of putative 
genetic clusters (K) and assign individuals into corresponding clusters. We conducted this analysis 
under admixture model and ran for 10 6 iterations with a burn-in length of 10 6 . The most likely K value 
was inferred by investigating mean Ln likelihood values. 

5. Conclusions 

In total, our study demonstrated that the cultured croaker had significantly reduced in genetic 
diversity in contrast to the wild and also presented statistically significant differentiation from the wild, 
indicating the current wild stock enhancement should be conducted with caution. With a relatively 
high level of genetic diversity, the wild populations showed important value for the ongoing breeding 
programs of this croaker. Simultaneously, we detected no differentiation among the wild populations, 
suggesting that the wild croaker in the northern and central East China Sea should be considered as one 
unit for management and conservation. 
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